Non-locality and short-range wetting phenomena 
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We propose a non-local interfacial model for 3D short-range wetting at planar and non-planar 
walls. The model is characterized by a binding potential functional depending only on the bulk 
Ornstein-Zernike correlation function, which arises from different classes of tube-like fluctuations 
that connect the interface and the substrate. The theory provides a physical explanation for the 
origin of the effective position-dependent stiffness and binding potential in approximate local theo- 
ries, and also obeys the necessary classical wedge covariance relationship between wetting and wedge 
filling. Renormalization group and computer simulation studies reveal the strong non-perturbative 
influence of non-locality at critical wetting, throwing light on long-standing theoretical problems 
regarding the order of the phase transition. 
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Density functional Q and interfacial Hamiltonian 
models 2] are complementary approaches to the theory 
of confined fluids. Mean-field, non-local density function- 
als give an accurate description of structural properties 
but are unable to account correctly for long-wavelength 
interfacial fluctuations. To understand these it is usually 
necessary to employ mesoscopic interfacial Hamiltonians 
based on a collective coordinate Z(x), measuring the local 
interfacial thickness. These models are essentially local in 
character containing a surface energy term proportional 
to the stiffness £ of the unbinding interface and a bind- 
ing potential function W(l). In more refined theories the 
stiffness also contains a position dependent term 3], 
which, it is has been argued, may drive the wetting tran- 
sition first-order JiJ. Despite progress over the last few 
years there are a number of outstanding problems par- 
ticularly for wetting with short-ranged forces. In addi- 
tion, recent studies of fluids in wedge- like geometries have 
uncovered hidden connections or wedge covariance rela- 
tions between observables at planar wetting and wedge 
filling transitions which have yet to understood at a 
deeper level. In this paper we argue that analogous to 
developments in density functional methods the general 
theory of short-ranged three-dimensional wetting should 
be formulated in terms of a non-local (NL) interfacial 
Hamiltonian. The model we propose directly allows for 
bulk-like correlations arising from tube-like fluctuations 
between the unbinding interface and the wall which 
contribute towards a binding-potential functional W[l, ip] 
where ip(x.) is the shape function for the wall. Whilst 
the possible importance of such tube- like modes has been 
muted by several authors || , to our knowledge this is the 
first quantitative treatment of them. Distinct contribu- 
tions to W[Z,^] reflect the number of interacting tubes 
and have a simple diagrammatic interpretation. Three 
implications of the NL model are discussed: A) In the 
limit of small interfacial fluctuations the NL model iden- 
tically recovers the known form of the local Hamiltonian 



allowing one to trace the specific position dependence of 
W(l), S(Z) to the Ornstein-Zernike (OZ) bulk correlation 
function. B) For filling in wedges with arbitrary tilt an- 
gles the model obeys the classical wedge covariance rela- 
tions observed in numerical studies of more microscopic 
theories. Such relations cannot be accounted for by lo- 
cal theories 0]. C) Renormalization group (RG) analysis 
show the non-perturbative influence of NL interactions. 
Despite precise connection with the model and RG the- 
ory of Fisher and Jin (FJ) at perturbative level the RG 
flow of the full NL model shows no stiffness instability 
and the wetting transition remains second-order. Simu- 
lations confirm these findings. 

Consider a Landau- Ginzburg- Wilson Hamiltonian 
based on a continuum order-parameter (magnetization) 
m(r) in a semi-infinite geometry with bounding surface 
described by a single- valued height function tp{x.) where 
x = (x, y) is the parallel displacement vector. Denoting 
the surface magnetization by TOi(x) we write 
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where ds^ = d\ + (Vip) 2 dx is the wall area element 
whilst 4>(m) and </>i(wi) are suitable bulk and surface 
potentials The bulk Hamiltonian is isotropic so the 
interfacial tension and stiffness are the same. Following 
FJ we identify the interfacial model H = Hlqw [to- (r)] 
where m=(r) is the profile which minimises Eq. (JQ sub- 
ject to a given interfacial configuration. FJ determined 
to h ( r ) perturbatively in terms of local planar constrained 
profiles Here we construct m~ non-perturbatively 
using Greens' functions or equivalently correlation func- 
tions defined within the constrained wetting layer. The 
latter reduces to the classical OZ form over relevant dis- 
tances provided the wetting layer is much thicker than 
its bulk correlation length k _1 . The NL Hamiltonian is 



H= /rfx|s v l + (V0 2 +Hl-Tp)\ +W[l,i>] (2) 
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which recovers the standard form of the binding potential 
appearing in local models. For the more general case of a 
non-planar interface near a planar wall, two contributions 
to the binding potential functional are local since 



nJ[Z,0] = J dsje-^ 1 ^, 



M =l,2 



(8) 



However fl 2 remains NL and can be rewritten as a two- 
body repulsive interaction 



FIG. 1: Schematic illustration of the diagrams which repre- 
sent the leading order contributions to W[l, rp]. 
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where h is proportional to the bulk field. There is no 
explicit position dependent tension but rather a binding 
potential functional with three leading contributions 

w[i,i/>] = -an 1 1 [i,i>} + b 1 n 2 1 [i,ip] + b 2 ttl[i,'i(>] (3) 

where a, b\ and b 2 are best regarded as phenomenologi- 
cal parameters to be identified later. Each fi^ represents 
integrated two-point interactions between /i and v points 
on the wall and interface mediated by the (rescaled) bulk 
OZ correlation function K(r) = Ke~ Kr /2iTr. These can 
be viewed as contributions to the free-energy of a con- 
strained thin-film arising from tube-like fluctuations of 
the bulk phase which tunnel from the interface to the 
wall 0]. The first term involves only one tube 
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etc and n% = 
is the distance between two 
The last two terms 
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involve two tubes and may be viewed as a self-interaction 
between points on the same interface or wall induced by 
the presence of a second surface. Each contribution can 
be represented diagrammatically as shown in Fig^ The 
upper and lower lines represent typical non-planar con- 
figurations of the interface and wall. The undulated line 
joining them represents the interaction function K (rxz) 
whilst the solid dots imply integration over the area of 
each surface. For general wall and interfacial configura- 
tions all contributions to W^,^] are NL. Simplifications 
arise when one or both are planar. If both the wall and 
interface are flat, ij){x) = 0, 1{x) = I the Hamiltonian per 
unit area W(l) — W[l,0]/A reduces to 



W(l) 



-ae 



+ b 2 e 



-2kI 



(7) 



where I = [Z(xi) + ^(x 2 )]/2 and S(x; I) is: 



S( X ;l) = - 
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valid for nl ^> 1. In the small gradient limit, |V/| <C 1, 
the NL term can be expanded and the model reduces to 



H[l,0] 



<lx i 5^ (vo 2 + W(l) 



with stiffness coefficient 



= S - aer Kl - 2b x nle 
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precisely recovering the FJ model and uniquely identify- 
ing a, b\ and b 2 - In particular a measures the deviation 
from the MF critical wetting temperature, b 2 oc a 2 and 
the sign of 61 determines the order of the MF transition. 
Thus the origin of the nle~ 2Kl contribution, crucial in 
the FJ analysis, can be traced directly to a perturbative 
treatment of the NL contribution ilf. 

Now consider fluid adsorption in a wedge geometry 
(-0 = tancv|a;|). The NL model satisfies the necessary 
requirement of classical wedge covariance known from 
numerical studies of the microscopic model (|TJ 6] . Clas- 
sical wedge covariance refers to the relationship between 
observables at MF critical wetting and MF wedge filling 
transitions. Let l n (0) denote the MF thickness of the 
critical wetting layer written as a function of the contact 
angle. Let l w (0, a) denote the thickness of the filling layer 
above the wedge mid-point. Numerical minimization of 
the LGW Hamiltonian shows that l w (9, a) = l w (0— a) 
as 9 — > a, for both shallow and acute wedges. This 
relation cannot be explained using an approximate lo- 
cal Hamiltonian in which the interface interacts with the 
(closest) wall via a binding potential dependent on the 
normal distance [10|. Such models predict the incorrect 
behavior l w (9,a) = secal. K (9 — a). In contrast the NL 
model obeys the correct wedge covariant relation. The 
reason for this can be traced to the structure of the NL 
binding potential. Since filling precedes wetting (a ^ 0), 
the dominant term is Now for a flat interfacial con- 
figuration Z(x) = Iq near a non-planar wall both fl\ and 



f2 2 are local with, for example 



l = ln 



J dx.y/1 + (V^) 2 e 



2„-K(io—0(x)) 
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showing that the effective local interaction occurs via the 
vertical distance to the surface. Near the filling phase 
boundary the interface is essentially flat in the filled sec- 
tion of the wedge and the fl\ contribution must be of the 
above form. This is sufficient to ensure covariance. We 
also remark that for wetting at more general non-planar 
walls the NL model reproduces the precise form of the 
stiffness matrix appearing in approximate two-field mod- 
els [HI valid for |VZ| < 1 and \Vip\ < 1. This means 
that in application to c omp lete wetting the NL theory 
satisfies exact sum rules |l2j . 

Finally we turn to the controversy surrounding fluctua- 
tion effects at planar critical wetting. The standard cap- 
illary wave (CW) model, obtained by setting E(Z) = E in 
Eq. f| 1 1 p . famously predicts non-universal criticality de- 
pendent on the wetting parameter lu — ksTn 2 /47rE |l3l . 
However, this strongly disagrees with fsing model simula- 
tion studies which show only minor deviations from 
MF-like critical wetting behavior (for the experiments, 
see Ref. 0])- The more refined F J model provides a pos- 
sible explanation of this discrepancy since the E(Z) term 
drives the transition first order for physical values of lu 
. Here we show that the stiffness instability is not a ro- 
bust mechanism since the wetting transition described by 
the NL model remains continuous. A linear RG theory 
can be constructed provided we first expand \Jl + (VZ) 2 
to square gradient order. The local terms fl^ generate ef- 
fective binding potential and position-dependent stiffness 
contributions which renormalize as in Refs. HEl. We 
focus on the renormalization of the NL potential S(x; I) 
which controls the order of the phase transition since it 
is responsible for the —le~ 2Kl term in the perturbative 
|VZ| <C 1 limit. After renormalizing up to a scale b = e* 
the NL term f2 2 retains its two-body form but with a 
modified potential St(x;l) satisfying the flow equation: 
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dt 
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l + J (Ax)\ c?St_ 
dt 
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where Jo (x) is a Bessel function of first kind and A is the 
momentum cutoff. This equation has the formal solution: 



S t (x; I) 



dl 



KSo(xe t ,l) exp ( 



-{kI-kI) 2 



4oj<I>(Aa;e t ,As) 



y / 4:TTUJ$(Axe t ,Ax) 



(15) 

where $(a, b) = ^ dt[l + J (t))/2t. We choose S (x, Z) = 
Q(l)S(x,l), with 0(Z) the Heaviside step function and 
S(x,l) given by Eq. (fTTTf) - As t — > oo, S t (x,l) becomes 
increasingly localized around x = 0. Using a match- 
ing technique we renormalize to a scale e at which the 
curvature of the effective binding potential Wt(l) at its 
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FIG. 2: (Color online). Plot of the mean wetting layer (Z) 
and surface magnetization operator Ami vs. h obtained by 
computer simulations of the CW, FJ and our NL model for 
lu = 0.8, a = 6 2 = and bi/K 2 k B T = 2.5. 

global minimum is of order of Ek 2 . Our Wt(l) has a local 
contribution due to the 0} and processes, and a NL 
contribution which is obtained from the expansion of f2 2 
in powers of VZ: 



W t NL (l) = 2tt / dxxS t {x,l) 



(16) 



Numerical integration of the RG flow equation show that 
the wetting transition is always second order, and quan- 
titatively similar to the non-universality exhibited by the 
CW model. This fact can be rationalized by noting that 
^(Axe 1 , Ax) ~ t as t — > 00 and kx < e _t , which is the 
range where So(xe t , Z) is non-negligible. Consequently, in 
our NL model there is no stiffness instability. The dif- 
ference with the RG predictions of the FJ model arises 
specifically from non-locality. Mathematically the FJ 
flow equations can be recovered from Eq. (|14l) if we ap- 
proximate the Bessel function term by its quadratic ex- 
pansion in x. However this is not valid at large distances 
and invalidates the stiffness instability. 

In order to check the RG predictions, we have per- 
formed Monte Carlo simulations of the CW, FJ and NL 
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Hamiltonians (with the approximation + (V?) 2 « 
1 + (V0 2 /2). Following Ref. 

we discretize by intro- 
ducing a Lx L lattice of spacing a with periodic bound- 
ary conditions in the directions parallel to the surface, 
but treating the interfacial position height as continu- 
ous variables. We chose a — 3.1623k _1 so that Ak _1 ~ 
7r / K<7 ~ 1) an d also set w = 0.8 and b\ — 2.5k 2 /cbT 
which are reasonable Ising-like parameters. We antic- 
ipate the critical wetting phase boundary remains MF 
(a = 0) for the CW and NL theories 0, whilst the 
FJ exhibits a first-order transition at higher tempera- 
tures 0| . Fig. [2] describes the behavior of the mean wet- 
ting layer thickness (I) and the surface magnetization-like 
operator Ami = (e~ Kl ) along the MF critical wetting 
isotherm a = 0, h — > 0. The FJ model clearly describes 
partial wetting in this limit consistent with a fluctuation- 
induced first-order transition. On the other hand the CW 
and NL models are qualitatively similar, showing contin- 
uous wetting. The divergence of the film thickness is 
well described by the RG result k{1) ~ — s/2ulnh even 
for moderately thick wetting layers. However the sur- 
face magnetization shows a much larger preasymptotic 
critical regime. The asymptotic non-universal behavior 
Ami ~ /i 1-1 / 2,y ii , with z/|| — ^/uj)~ 2 is not observed 

until the wetting layer k{1) ~ 10 for very large lattice 
sizes kL ~ 300. This is strongly suggesting that current 
Ising model simulations will not be able to observe sig- 
nificant deviation from MF behavior provided they focus 
on surface quantities. 

We conclude by mentioning extensions and limitations 
of our theory. It is straightforward to generalize the 
model to describe wetting at heterogeneous substrates 
with, for example, hydrophilic and hydrophobic domains 
[l7| . For such geometries all the contributions to W[l, ip] 
will be NL and may influence mesoscopic wetting behav- 
ior. The same is true for other types of homogeneous 
sculpted substrates as, for example, wedges of parabolic 
cross section [l8j. However, even for the simple case 
of wetting at homogeneous planar walls, systematic im- 
provements of the theory can be envisaged. The binding 
potential functional @ is only valid for wetting layers 
many times larger than kT 1 . If the wetting layer is only 
a few bulk correlation lengths thick, the interaction func- 
tion K (r) should be replaced by the correlation function 
defined within the constrained profile. This would mod- 
ify the form of the binding potential functional at short 
distances and may influence how the asymptotic critical 
regime is approached. 

A.O.P. gratefully acknowledges Prof. D. B. Abraham 
for discusssions regarding the origin of the binding poten- 



tial. We also thank Prof. M. E. Fisher and Dr. J. Gib- 
bons for helpful comments. J.M.R.-E. acknowledges par- 
tial financial support from Secretaria de Estado de Ed- 
ucation y Universidades (Spain), co-financed by the Eu- 
ropean Social Fund, and from the European Commission 
under Contract MEIF-CT-2003-501042. A.L. acknowl- 
edges partial financial support from the A. G. Leventis 
Foundation. 



[1] See for example, R. Evans in Fundamentals of Inhomo- 
geneous Fluids, edited by D. Henderson (Dekker, New 
York, 1992). 

[2] G. Forgacs, R. Lipowsky and Th. M. Nieuwenhuizen, in 
Phase Transition and Critical Phenomena, vol 14, edited 
by C. Domb and J. L. Lebowitz (Academic Press, New 
York, 1991). 

[3] M. E. Fisher and A. J. Jin, Phys. Rev. B 44, 1430 (1991); 

A. J. Jin and M. E. Fisher, Phys. Rev. B 47, 7365 (1993). 
[4] M. E. Fisher and A. J. Jin, Phys. Rev. Lett. 69, 792 

(1992); A. J. Jin and M. E. Fisher, Phys. Rev. B 48, 

2642 (1993); C. J. Boulter, Phys. Rev. Lett. 79, 1897 

(1997). 

[5] A. O. Parry, M. J. Greenall and A. J. Wood, J. Phys.: 
Condens. Matter 14, 1169 (2002). 

[6] M. J. Greenall, A. O. Parry and J. M. Romero-Enrique, 
J. Phys.: Condens. Matter 16, 2515 (2004). 

[7] D. B. Abraham, J. T. Chayes and L. Chayes, Commun. 
Math. Phys. 96, 439 (1984). 

[8] D. B. Abraham, private communication. See also M. E. 
Fisher, in Statistical Mechanics of Membranes and Sur- 
face, edited by D. R. Nelson, T. Piran and S. Weinberg 
(World Scientific, Singapore, 1989). 

[9] H. Nakanishi and M. E. Fisher, Phys. Rev. Lett. 49, 1565 
(1982). 

[10] K. Rejmer, S. Dietrich and M. Napiorkowski, Phys. Rev. 

E 60, 4027 (1999). 
[11] C. J. Boulter and A. O. Parry, Phys. Rev. Lett. 74, 3403 

(1995). 

[12] See J. R. Henderson in Fundamentals of Inhomogeneous 
Fluids, edited by D. Henderson (Dekker, New York, 
1992). 

[13] E. Brezin, B. I. Halperin and S. Leibler, Phys. Rev. Lett. 

50, 1387 (1983); D. S. Fisher and D. A. Huse, Phys. Rev. 

B 32, 247 (1985). 
[14] K. Binder, D. P. Landau and D. M. Kroll, Phys. Rev. 

Lett. 56, 2272 (1986). 
[15] D. Ross, D. Bonn and J. Meunier, Nature 400, 737 

(1999). 

[16] G. Gompper and D. M. Kroll, Phys. Rev. B 37, 3821 
(1988). 

[17] H. Gau et al, Science 283, 46 (1999). 

[18] C. Rascon and A. O. Parry, Nature 407, 986 (2000). 



